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We present the results of our perturbative calculations of the static quark potential, small Wilson loops, the 
static quark self energy, and the mean link in Landau gauge. These calculations are done for the one loop 
Symanzik improved gluon action, and the improved staggered quark action. 



It has long been known that perturbative calcu- 
lations are necessary for precision measurements 
of many quantities. Unfortunately perturbation 
theory using a lattice cutoff is difficult, due to ex- 
tra diagrams and vastly more complicated Feyn- 
man rules. The complexity of these calculations 
is an impediment to doing the higher-order per- 
turbation theory required for many important ap- 
plications of improved actions. The present work 
attempts to automate as much of the perturba- 
tion theory as possible in order to make these 
types of computations more straightforward. 

For any given action some of the most basic 
things we would like to compute are various op- 
erators that can be built out of products of links. 
These include the static quark potential, small 
Wilson loops, the static quark self energy, and 
the mean link in Landau gauge. 

Perturbative determinations of these quantities 
are central to many other applications. For ex- 
ample, the static quark potential can be used to 
determine a physical, renormalized coupling, and 
the perturbative expansions of small Wilson loops 
can be used to extract the strong coupling from 
simulations (see [Q). We have the techniques in 
place to calculate all of these, for an arbitrary 
gluon action, along with the contributions from 
quark loops for highly improved quark actions. 
This paper reports preliminary results for the 
above quantities, for the one-loop Symanzik im- 
proved gluon action, and the improved staggered 
quark action. 



These calculations are similar to traditional 
continuum perturbative QCD. The major differ- 
ence is that the ultraviolet regulator is the lat- 
tice cutoff, which leads to the computational dif- 
ficulties mentioned above. A further problem is 
how to regulate infrared divergences beyond one- 
loop, in a gauge covariant manner. This can be 
done using twisted periodic boundary conditions 
1^. Feynman rules and factors from the opera- 
tors are generated using the Liischer and Weisz 
vertex generation algorithm PYTHON or 

CH — h scripts implement this program for a given 
set of links. An advantage to this method is the 
ease with which new actions and operators can be 
treated. We use FORTRAN programs to do the 
momentum sums. These are generally done on 
large but finite lattices (typically 100'* volume). 

The action we consider here is the one-loop 
Symanzik improved action for isotropic lattices, 

Sg = I3pi 51 (1 - ^M'^) + Prt 51 (1 - ^M-) 
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where 

Ppi = ^, /3rt = -T^(l + 0.4805a«) 

/3p„ = -^0.03325as- (2) 

ul 

We present our results for an expansion in the 
bare lattice coupling aiatt — 5^/(47r). We include 
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tadpole counterterms for (|l|), which are generated 
by the first order expansion of the mean-field in 
the rectangle terms. In this paper we use the 
mean-field, uq, defined by the average plaquette 
^ which, for the action of (|l]), is given to first 
order by uq ~ f — 0.7671aiatt- Our second order 
results also included the counterterm generated 
by the 1x1x1 paths in (|l]). 

Next we consider the static quark potential 
which is the central quantity that needs to be 
computed. We can use it to define a renormalized 
coupling ay, which can be used as the expansion 
parameter for other quantities. 

To compute the static quark potential we take 
the correlator of two Polyakov lines of length Lt, 
separated by a distance R, < L{R,Lt) >■ The 
static quark potential is then given by 

V{R) + 2Eo= lim ^ \n < L{R, Lt) > , (3) 

where Eq is the self energy of an isolated quark 
(see below for our calculation of this quantity). 
Expanded in the bare coupling the static quark 
potential should have the following form [Q 

x{l + ax^tt{2l3oln{T:R) + C{R)), (4) 

where /3o = (H - 2/3n/)/(47r). 

As mentioned above, the static quark potential 
is used to define a renormalized coupling. This is 
done by demanding that the Fourier transform of 
(0) have the form 



V{q) 
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Using this definition, we obtain the expansion for 
the bare coupling in terms of the physical one. 
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As a test of our calculations we reproduced 
the known result for the Wilson gluon action, 
C = 4.702. We have also determined C for the 
Symanzik improved gluon action. We find 



We could easily recompute C for other gluon (and 
quark) actions, and with other definitions of uq. 

In addition to the static quark potential, we 
have computed a number of other quantities. Ta- 
bles and H give results for the logarithms of 
small Wilson loops, whose perturbative expan- 
sion is defined by, 
1 
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w^{R,T)aU. (8) 



The results for the Wilson action agree with those 
of §. 

We have also computed the static quark self 
energy Eq through second order. We define the 
self energy Eo{L) on a finite lattice according to 



ML) = -^HPtiL)] = ^c„(iXtt. 



(9) 



Ci = 3.23(13). 



(7) 



Here Pt (L) is the Polyakov line on a lattice of size 

Figures |^ and ^ show our results for the first 
and second order coefficients on a series of vol- 
umes, for three sets of boundary conditions (Wil- 
son glue, PBC stands for periodic boundary con- 
ditions, Txy and Txyz stand for twisted peri- 
odic boundary conditions along two and three 
spatial planes, respectively). The infinite vol- 
ume extrapolation agree with earlier estimates, 
Eq = 2.1172aiatt + 11.152ai2^jj. These finite vol- 
ume results were used in a determination of the 
third order self energy, using Monte-Carlo meth- 
ods Q . We have also calculated the self energy 
for the improved gluon action (m) with the result: 
£^0 = 1.8347(5)aiatt + 8.01(2)^^^^. 

Finally, we have computed the mean link in 
Landau gauge. In agreement with earlier deter- 
minations we have for the Wilson action, uq = 
1 - 0.9738(2)aiatt - 3.33(l)Qf^jj. For the one loop 
Symanzik improved action we report, uq = 1 ~ 
0.7501(l)aiatt-2.06(l)a2^tt. 

The methods for automatic vertex genera- 
tion can also be readily applied to complicated 
fermionic actions. For example, we computed the 
Uf part of the average plaquette at second or- 
der for improved staggered fermions . We find: 
W2(l, 1) = 1.958(2) - 0.06969(4)n/. Calculations 
of the second order n f pieces for the other quan- 
tities considered in this paper are in progress. 
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Figure 1. The C'(aiatt) contribution to the static 
quark self energy. 
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Figure 2. The 0{ai^^^) contribution to the static 
quark self energy. 
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We turn now to some conclusions. Our cur- 
rent results demonstrate the versatility of our 
approach. High precision results, necessary for 
accurate determinations of many quantities, are 
currently being generated. For improved stag- 
gered quarks this includes one loop improvement 
of the action, and the matching of the quark cur- 
rents and four quark operators |^ as well as quark 
mass renormalization Bl. 



Table 1 

Perturbative Wilson loops evaluated using Wil- 
son glue, errors are from the VEGAS integrations. 
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Table 2 

Perturbative Wilson loops 
Symanzik improved glue. 



evaluated using 
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